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The dominance of S and D pairs in the description of deformed nuclei is one of the facts that 
provided sustain to the Interacting Boson Approximation. In Ref.pJ, using an exactly solvable 
model with a repulsive pairing interaction between bosons it has been shown that the ground state 
is described almost completely in terms of S and D bosons. In the present paper we study the 
excited states obtained within this exactly solvable hamiltonian and show that in order to obtain a 
rotational spectra all the other degrees of freedom are needed. 
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INTRODUCTION 

The interacting boson approximation (IBA)0,0,Q as- 
sumes that the elementary excitation used for the de- 
scription of the low energy nuclear spectra in open shell 
nuclei are formed by pairs of nucleons coupled to an- 
gular momentum and 2 (S and D). These excitations 
behave in this approximation as bosons. An intriguing 
question is how to determine the structure of these S and 
D bosons: 

It seems natural to use the Hartree-Fock-Bogoliubov 
(HFB) approximation where the ground state wave func- 
tion can be written as in Ref. Q. 



|$ >= Nexpi-J^Zkiala}), 



(1) 
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where at creates a particle in a state k, N is a normal- 
ization constant and Z is the antisymmetric structure 
matrix 0, 0- 

In order to study this point it may be enough to use 
the results obtained with schematic calculations using 
the pairing plus quadrupole hamiltonian in the Nils- 
son+BCS approximation. These calculations show that 
the S and D bosons are dominant in the description of 
the ground state if it is considered a coherent state of 
Cooper pairs 8] or as a pair condensate 0. However, self 
consistent HFB calculations for the sing le j-shell[irj or 
many shells 11, 12] using unrenormalized simple interac- 
tions show that, in these cases, the space spanned by the 
S and D pairs does not provide a reasonable description 
of well-deformed systems. In Ref. 01 it is shown that if 
the interaction strengths are renormalized to reproduce 
the value of (3 (deformation parameter) and A (the super- 
conductive gap) a reasonable description of the ground 
state of wcll-dcformcd nuclei is obtained. 



In Ref.flJ] a large group of integrable many body 
hamiltonians (as well as the constants of motion related 
to them) was found. One of these hamiltonians is a 
repulsive pairing interaction between bosons with single 
boson energies that may have any value. For a particular 
set of single particle energies and a repulsive pairing 
interaction they found that for a large number of bosons 
(N) the ground state of the system only contains two 
types of bosons (in that paper the relevant bosons were 
S and P, all other bosons have occupations numbers of 
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The reason why it is interesting to consider 



this hamiltonian is that the repulsive pairing interaction 
between the bosons mucks up the influence of the Pauli 
principle between the quasibosons built on fermion pair 
(see Appendix 1 for a description of this relation in the 
single j-shell case). In a recent letter Q a renormalized 
pairing interaction as well as a set of single bosons 
energies more related to the nuclear case yields that 
the S and D states are the only ones that have large 
occupations numbers. In Ref.|lq] a similar calculation 
is done using group theoretical methods (therefore they 
are forced to use the pairing interaction instead of the 
renormalized one) and the S-D dominance is attributed 
to the choice of boson's energies. 

The present paper uses the recently proposed 
Richardson-Gaudin (RG) integrable models 14] 
(for a recent review see |16() to study the excitations 
obtained when using a repulsive pairing interaction 
between bosons for "realistic" situations. By "realistic" 
situations we mean the use of single boson's energies 
similar to the ones that correlated pairs of nucleons 
have in nuclei. The strength of the pairing interaction 
is used as a free parameter. It must be noted that for 
the single j-shell the value of the pairing interaction 
between the "bosons" is minus the pairing interaction 
between fermions, as shown in Appendix 1. The use of 
this treatment allows for the study of two problems: the 



2 



first one is the influence of using realistic values for the 
boson energies on the S-D dominance as well as for the 
interaction strength; the second one is the possibility of 
constructing a rotational spectra if a multiboson space 
is considered. 



HAMILTONIAN USED AND ITS TREATMENT 

The hamiltonian that will be used is 

H = E + f £(^r + k^k+), ( 2 ) 



where ei is the energy of the boson I and the operators 
= and K~ are the generators of an SU(1,1) algebra. 
The set of generators Kf = \hi + Kf and K~ = 
(K^Y satisfies the commutation relations 



[Kt,Kj]=-25 l3 Kl 



(3) 



This SU(1,1) algebra is realized in terms of operators 
that create and annihilate bosons. In the usual pseudo- 
spin representation, the generators can be written as 



j 2 / J 3 



j i 1 ^ — ^ 

— > Kj — 2 a jm a jm 



,t J 



(4) 



where aj m (a jm ) creates (annihilates) a boson in the state 
\jm), \jfh) is the state obtained by acting with the time 
reversal operator on the state \jm), and fl, is the total 
degeneracy of the single-boson level j. Here the realiza- 
tion involves correlated pair operators. It is simple to 
express the pairing hamiltonian in terms of these opera- 
tors. 

The Hilbert space of N particles moving in shells 
k,l,m,....p can be classified in terms of the product 

of groups SU(l,l)fcxSU(l,l) ,xSU(l,l) m x SU(1,1)„. 

All the states of the Hilbert space can be written as 



., n p v >- 



1 



(K+rK...(K+r»\v >, (5) 



where Af is a normalization constant. The possible num- 
bers of pairs in each level is < n s < N and the states 
\v >= |z/fe,i/;, v m , ...v p > of unpaired particles are defined 

as 



There are three families of fully integrable and exactly- 
solvable RG models that derive from these SU(1,1) alge- 
bras, which are referred to as the rational, trigonometric 
and hyperbolic models, respectively 

Hi- 

The pairing 

interaction is related to the "rational" family, for which 
they express the integrals of motion in terms of the above 
generators as 



Ri = Kf + 2gJ2 ^ —L— [KfKr+Kr K f} 



»'(^J) 



K°K° 



(7) 



For each degree of freedom i, there is one real arbitrary 
parameter that in the rational model is equal to e, that 
enters the integrals of motion. It can be easily checked 
that these operators commute among themselves, and 
each one commutes with the operator K° = K?- 

The eigenvalue equation = r^i) can be readily 

solved using methods analogous to those first introduced 
by Richardson to treat the quantum pairing problem |l7| . 
The eigenvectors take the form 



M 



7 T K l> 

(e a - ei) 



(8) 



where \v) is a state that is annihilated by all the K~ 
and M is equal to the number of collective B + operators 
that comprise the state. The structures of the collective 
operators are determined by a set of M parameters e a , 
which satisfy the set of coupled nonlinear equations 



l + 4 5 ^-^+4 5 ]T ^^=0. (9) 



The associated eigenvalues take the form 



{ J a ) 



(10) 



We note here that each independent solution of the set 
on nonlinear coupled equations @ defines an eigenstate 
(JSJl and eigenvalues that are given by 



(11) 



K-\v>=0, 



nJu >= vAv >, 



K°\v >= d„ 



\v> 
(6) 

where d s = \{2v s + Q s ) and N = 2M + v, N being the 
total number of bosons, M the total number of pairs and 
v (usually called the total seniority) the total number of 
unpaired bosons that is the sum of all the v s . 



Following Richardson |17| , the occupation numbers can 
be calculated as 



8H P \ _ ^ de a 



(12) 
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From (JSJ, and making use of © and (fTT|l . a set of M 
coupled linear equations in terms of M new unknowns 
can be obtained, which give the A boson occupation num- 
bers. Details of the derivation can be found in ref. [l9j . 
Dividing these occupation numbers by the total number 
of bosons we obtain the occupation probabilities of the 
different bosons. 



RESULTS AND DISCUSSION 

In order to relate the calculation to the nuclear case 
it is necessary to choose the properties of the boson de- 
grees of freedom for each angular momentum as well as 
the coupling constant. The microscopic image that we 
have in mind is the Nilsson model and therefore it seems 
natural to consider only one type of bosons for each even 
angular momentum. For angular momentum two the se- 
niorities needed to describe the low lying states are well 
known from the study of the five dimensional har- 
monic oscillator. For A = 2 we will consider only states 
with seniority one (3=2), two (J=2 and 4) and three 
(J=0, 3,4,6)). A similar treatment can be used for the 
A bosons (the first excited states corresponding to these 
degrees of freedom will have seniority one and angular 
momentum J = A and are the only ones that we will 
consider). It is worthwhile to remark that as the Hamil- 
tonian has vanishing matrix elements between states with 
different seniorities, the states connected by the hamil- 
tonian must have the same seniorities. There are too 
many parameters if the boson's single particle energies 
are left as free parameters. We performed all our cal- 
culation with a large but meaningful number of bosons 
(40) and we choose the single boson energies by imposing 
the condition that the ratios of calculated energies of the 
first 2+, 4 + , 6 + and 8 + states are as close as possible to 
the values provided by the rigid rotor for G/hv2 = .72 
(for this value of Gfhv^ the ratios of the energies of the 
ground state band have reached what we can consider 
their asymptotic values) taking into account bosons with 
all the even angular momentum up to 12. As we do not 
know the ground state energy, we can only determine the 
differences between the bosons single particle energies. 
We define our energy scale by making the difference of 
energy between the A = 2 and A = bosons equal to 
one. The resulting energy differences (to the J=0 boson) 
have values which are reasonables: E4 = 1.25; Ee = 1.50; 
E 8 = 2.12; Em = 2.28 and ; E X2 = 2.42. 

In Fig. 1 we show the excitation energy of the 2i state 
as function of G/hv2 (being hv2 the energy of the A = 2 
boson). The asymptotic value (w .124) provides us with 
an absolute energy scale. If this state is the 2i state of 
the ground state rotational band the energies are of the 
right order of magnitude and G/I1V2 must have a value 
of the order of 0.2 to reproduce the fcrmionic strength. 

Initially we include all bosons with even angular mo- 
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FIG. 1: Energy of the 2i state as function of the scaled 
strength parameter G/hv2- 
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FIG. 2: Ratio of the energy of the Ji states to the energy 
of the 2i state as function of the scaled strength parameter 
G/hv2 for a calculation that includes all even values of J up 
to 12. The rigid rotor values are shown to guide the eye. 

mentums up to A = 12. In Fig. 2 we show the ratio 
between the energy of the first state for each even angu- 
lar momentum (up to J=12) to the 2i state as function 
of G/hv%. It is evident the departure from the rigid ro- 
tor values when the angular momentum increases or G 
decreases. In the region of Gfhvi « 0.2 the structure 
of the excitations is similar to the rigid rotor. For the 
model we are using, the states are labeled by the senior- 
ity of the different bosons. To get a feeling of the types 
of result obtained in table 1 we list the excitation ener- 
gies for G/hv2 = 0.72 for some selected sates, showing 
also the possible angular momentum of these states. It 
follows that there are many excited levels in the low en- 
ergy region and that it is difficult to recognize the excited 
rotational bands. 
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FIG. 3: Occupation probabilities of the different A bosons as 
function of the scaled strength parameter G/hi/2 for a calcu- 
lation that includes all even values of A up to f 2. The lines 
corresponding to sizeable contributions are labelled by A. 



In Fig. 3 we show the occupation probabilities of the 
different boson levels obtained for some selected states. 
In general only two boson levels are populated (the S 
and D bosons share more than 85% of the occupation 
probabilities). The contribution of the A = boson will 
be more important if we increase the number of bosons 
or if we consider a renormalizcd interaction as the one 
used in Ref.0. 

In Table 2 we show the occupation probabilities for the 
different boson degrees of freedom for G/hv2 = 0.72 and 
for the first state corresponding to different seniorities. 
From this table follows that the contributions of the 
different bosons (with seniority one) are important in 
the description of the ground state band, even if the 
occupation probabilities are small. 

A question that one may ask is how will change the de- 
scription obtained as we diminished the number of boson 
degrees of freedoms. For that purpose we also performed 
the calculation including only bosons up to A = 8. We 
tried again to reproduce as well as possible the rigid rotor 
by changing the bosons energies. The values obtained for 
the energies (in the region of large G) were E4 = 1.10; 
E6 = 1.37; Eg = 2.11 which are similar to the values used 




FIG. 4: Ratio of the energies of the J 1 states to the energy 
of the 2i state as function of the scaled strength parameter 
G/hv2 for a calculation that includes all even values of A up 
to 8. The rigid rotor values are shown to guide the eye. 



in the larger calculation. The asymptotic value for the 
2 + energy was 0.101. 

In Fig. 4 we show results similar to Fig. 2 but using 
the smaller space. It is clear that we are not able to 
reproduce properly the 8 + state which gives an indication 
of the importance of the A = 10 and A = 12 bosons in 
its description. In Figure 5 we display the occupation 
probabilities for this case and in Table 3 we reproduce 
the equivalent of Table 2. 

In Fig. 6 we show the energy ratios obtained for the 
ground state band when we consider only S and D bosons. 
In this case there are not free parameters (except for the 
pairing strength) and it is quite clear that we are not able 
to reproduce a rotational band. 

The main weakness for the relevance of the present cal- 
culation is related to how well is taken into account the 
Pauli Principle by our model hamiltonian. In a Nilsson 
like calculation one has two main ingredients: on the one 
side the quadrupole-quadrupole interaction deforms the 
systems, whose ground sate can be written as a conden- 
sate of a boson that if formed by S, D, G, etc. bosons 
with well-defined probabilities 0, and then one must 
consider the pairing interaction, that change by a small 
amount these probabilities. The contribution due to the 
pairing interaction between fermions can be taken into 
account using Pauli violating diagrams 0. The dia- 
grams that appear in this treatment are similar to the 
ones generated by our model hamiltonian, but they do 
not have constant matrix elements. The fact that one 
obtains a rotational spectra, when using a large space, is 
an indication that the replacement of the correct matrix 
element by a constant may not be a bad approximation. 

We can now go on to the main conclusion of the present 
paper. The S and D bosons used in the IBA model 0iQ 
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FIG. 5: Occupation probabilities of the different A bosons as 
function of the scaled strength parameter G/hv2 for a calcu- 
lation that includes all even values of A up to 8. The lines are 
labelled by A. 



are different from the ones that appears in a quadrupole- 
quadrupole+pairing (Nilsson like) calculation. Using our 
model hamiltonian we have shown that one cannot ignore 
the importance of the high angular momentum bosons in 
the description of the ground sate rotational band, even 
if the occupation probabilities of these bosons are very 
small. Our result suggests that the S and D bosons used 
in the IBA are not directly related to the ones obtained 
in a Nilsson+BCS treatment. The S and D dominance 
obtained is a manifestation of a general effect found in 
Ref. , but the small amplitudes of the "bosons" with 
large angular momentum are essential to obtain a good 
description of the ground state rotational band. 
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APPENDIX 1 

The single j shell 
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FIG. 6: Ratio of the energies of the states with seniority v to 
the energy of the 2i state (corresponding to seniority one) as 
function of the scaled strength parameter G/hv2 for a calcu- 
lation that includes only S and D bosons. 



An attractive pairing-interaction between fermions in a 
single j-shell has a well known exact solution that for the 
seniority zero states is expressed in terms of the number 
of pairs of particles (N) and the degeneracy (Q = j + |) 
as 



E(N) = — \G\Nft 1 



(N-l) 



(13) 



The second term can be considered to arise from the 
Pauli principle between the pairs of fermions, each having 
an energy — |G|f2. 

In the case of S bosons the repulsive pairing hamilto- 
nian is written as (G > 0) 



Hp = WsS^S + GS' l S jl SS, 



(14) 



and it is diagonal in the states of N bosons, yielding for 
the energy 

G(N - 1) 



E(N) = W S N + GN(N - 1) = W s N 1 + 



(15) 

To relate this result with the fermionic case one must 
use W s — — \G\fl and then both energies coincide. 
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TABLE I: Excitation energies, in units of hv2, for the three first solutions with the seniorities displayed in the table. A list of 
all the angular momentum related with this seniority structure are given in the last column. 



M 


^J=0 


V.J=2 








E 1 


E'2 E$ 


Possible angular momentums 


20 

















0.00000 


0.59667 1.35124 





f9 


1 


1 











0.12424 


0.79629 1.51183 


2 


f9 


1 





1 








0.41361 


0.99778 1.79228 


4 


f9 


1 








1 





0.86824 


1.46588 2.18733 


6 


f9 


1 











1 


1.54591 


2.14410 2.88702 


8 


f9 





2 











0.29651 


1.03966 1.72393 


2,4 


18 


1 


3 











0.51105 


1.32094 1.98002 


0,3,4,6 


18 


1 


2 


1 








0.75453 


1.50430 2.21865 


0, 1,2 2 ,3 2 ,4 2 ,5 2 ,6 2 ,7,8 


18 


1 


2 





1 





1.19265 


1.95006 2.59798 


2,3,4 2 ,5 2 ,6 2 ,7 2 ,8 2 ,9, 10 


18 





3 


1 








0.993228 1.81334 2.50029 0, l 2 , 2 3 , 3 3 , 4 4 , 5 3 , 6 3 , 7 3 , 8 2 , 9, 10 



TABLE II: Occupation probabilities for the bosons with angular momentum J of the first states of a given seniority for the 
calculation that includes all even J up to 12 with G = 0.72. In the last column are indicated the structure and angular 
momentums of the corresponding states 

P(J=0) P(J=2) P(J=4) P(J=6) P(J=8) P(J=10) P(J=12) States 



0.70384 0.17515 


0.06277 0.02505 


0.01092 


0.01104 


0.01124 


J=0 {vi = 


0) 


0.69740 0.19038 


0.05658 0.02372 


0.01048 


0.01062 


0.01082 


J=2 (i/2= 


1) 


0.69807 0.14543 


0.09903 0.02464 


0.01079 


0.01092 


0.01112 


J=4 (>4 = 


1) 


0.69602 0.15579 


0.06154 0.05365 


0.01085 


0.01098 


0.01118 


J=6 (y 6 = 


1) 


0.69431 0.15963 


0.06183 0.02489 


0.03715 


0.01099 


0.01120 


J=8 (v 8 = 


1) 


0.69405 0.16008 


0.06186 0.02489 


0.01087 


0.03704 


0.01120 


J=10(z/i = 


= 1) 


0.69386 0.16040 0.06189 0.02490 0.01087 


0.01100 


0.03710 


J=12(l/12 = 


= 1) 


0.68952 0.20713 


0.05060 0.02228 


0.00999 


0.01014 


0.01034 


J=2,4 (iaj= 


=2) 


0.68055 0.22453 


0.04514 0.02082 


0.00964 


0.00964 


0.00984 


J=0,3,4,6 (iaj = 3) 
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TABLE III: Occupation probabilities for the bosons with an- 
gular momentum J of the first states of a given seniority for 
the calculation that includes all even J up to 8 with G = 0.72. 
In the last column are shown the structure and angular mo- 
mentums of the corresponding states 

P(J=0) P(J=2) P(J=4) P(J=6) P(J=8) States 

0.64192 0.22291 0.08728 0.03680 0.01108 J=0 (yi = 0) 

0.63763 0.23677 0.07959 0.03525 0.01078 J=2 (i*2=l) 

0.63747 0.18653 0.12862 0.03638 0.01100 J=4 {u 4 =l) 

0.63523 0.20049 0.08600 0.06723 0.01105 J=6 (v 6 =l) 

0.63270 0.20684 0.08643 0.03667 0.03736 J=8 (z/«=l) 

0.63212 0.25220 0.07181 0.03346 0.01041 J=2,4 {u 2 =2) 

0.62563 0.26826 0.06451 0.03158 0.01002 J=0,3,4,6 (u 2 =3) 



